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>  What  is  the  Computer  Simulation  (Modeling) 

>  Why  Do  We  Need  Simulation 

>  Computer  Modeling  -  Basics: 

4-  Experimental  Characterization 
4-  Process  Modeling 
4  Analytic 
4  Numeric 

>  Computer  Modeling  -  Advanced: 

4-  Control  and  Sensing 
4-  Optimization 

>  Putting  the  Whole  Package  Together  or  the  Art  of 
Integration 

>  Conclusions 


Simacek  ONR  Workshop  -  3 


©2003  University  of  Delaware  All  rights  reserved 


2  July  2003 


RTM  and  VARTM  Processes 


RTM 


5.  Demolding  and  Final 
Processing 


1.  Preform  Manufacturing 

o 


2.  Lay-up  and  Draping 


4.  Resin  Injection  and  Cure 


VARTM 


1.  Preform  Manufacturing 

o 


2.  Lay-up  and  Draping 

o 


y  y  y  y 

o 


3.  Bagging 


AAA  A 


6.  Demolding  and  Final 
Processing 


5.  Resin  Injection  and  Cure 


o 


4.  Drawing  Vacuum 
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What  Is  To  Be  Simulated? 
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The  Model 


Input(s) 

Material  Parameters 

•  Geometry 

•  Resin  Viscosity 

•  Preform  Porosity 

•  Permeability 

•  Environment 

•  Global  Temperature 


Process 


Output 


Governing  Equation(s) 


Numerical  Solution 


Process  parameters 

•  Time 

•  Pressure 
Part  Parameters 

•  Flow  Patterns 

•  Saturation 
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Model  Creation 


>Choose  the  output  of 
interest 

4-  Decide  the  accuracy 

> Figure  out  the  governing 
equations 

4-  Decide  of  what  is 
important 

4-  Decide  of  what  input  and 
environment  parameters 
are  considered 

> Solve  the  equations 

4-  Numerical  solution 

>  Number  of  Sources  of  Error 

4-  Mathematicians  mislead! 
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Outpm 


Material  Parameters 
♦Geometry 
♦Resin  Viscosity 
♦Preform  Porosity 
♦Permeability 
Environment 
•Global  Temperature 


Governing  1 

Equation(s) 

Numericc 

il  Solution 

Process  parameters 

♦Time 

•Pressure 

Part  Parameters 

♦Flow  Patterns 

♦Saturation 


Keel  Beam  And  Tooling 
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Flow  Fronts  in  a  40-ft.  Keel  Beam 
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Flow  Fronts  in  a  40  ft  Keel  Beam 
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Why  Do  We  Need  Simulation  Tools? 


Can  We  See  Into  the  Mold? 

Can  We  Measure  What  is  Happening? 

If  We  Could,  How  Much  Good  Can  Come  Out  Of  It? 
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Can  We  See  The  Flowfront? 


Stiffening  bar 


Simacek  ONR  Workshoi 


Can  We  ‘Measure’  it? 


Sensor 

wires 


Flexible 

Tekscan 

sensor 


PSI 

255.00 

198.46 

179.62 

160.77 

141.92 

125.08 

104.25 

85.58 


M-nn  hH 
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Force  sensitive  layer 
(changes  resistance) 

Cover  layers 
(MYLAR) 


Handle 
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In  Order  to  Obtain  Visual  Input,  You 
Must  Have: 
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In  Order  to  Obtain  Input  From 
Simulation,  You  Need: 
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Why  Do  We  Need  Simulation  Tools? 


>  Can  We  See  Into  the  Mold? 

4-  In  RTM  Manufacturing  NO! 

4  In  VARTM  or  Laboratory  A  little 

>  Can  We  Measure  What  is  Happening? 

4-  Yes,  but  at  High  Cost.  On  Top,  the  Used  Technology  is 
Immature. 

>  If  We  Could,  How  Much  Good  Can  Come  Out  Of  It? 

4-  Visual/Sensorial  Examination  Requires  Mold  and  Part  to  Be 
Real 

4-  Errors  ->  Post  Mortem  Analysis  ->  New  Trial 

>  PREDICTIVE  CAPABILITY  IS  NEEDED! 
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What  Is  Needed  for  Modeling 


>  Computer  Modeling  Basics: 

♦  Experimental  Characterization  of  Used  Materials. 

♦  Process  Modeling  -  Equations  and  Their  Solution. 

4  Analytic. 

4  Numeric. 

4-  Tricks  of  the  Trade. 

>  Computer  Modeling  Advanced: 

4-  Modeling  Control  and  Sensing  and  “Intelligent” 
Manufacturing  Processes. 

♦  Process  Optimization. 
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Lesson  I. :  Knowing  the  Material 


>  Characterization  of  Materials 

4-  Resins 
4-  Preforms 

>  Preform  Data  that  Influence  the  Process 
4-  Permeability 

4-  Porosity 

4-  Thermal  Characteristics 

>  Ways  to  Obtain  Them 

4-  Computer  Simulations  (Hypothetic) 

4-  Experimental 

>  Other  Uses  of  Experiment 

4-  Verification! 
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Numerical  Permeability  Prediction 


z  A 


I 


Mold  wall  or  another  layer 


Resin  channel 


►  X 


Warp  tow 

l-D  (in  z  direction)  Darcian  flow 

Weft  tow 

In-plane  (x-y),  lubrication  flow 


K„(pB  -  pA)  1  /  dUhWdgfi/dx)  c)((hA)3d  pfyd  v)\  n 

nhF  6r|  \  dx  *  dy  J~U 

K77(pB -  pA)  1  /  9((hB)3dpB/dx)  8((hB)39pB^y)\ _  n 

TihF  “  6ri  \  dx  +  dy  )-" 


Input  to  model:  geometry:  hA(x,y),hF(x,y),hB(x,y )  Output :  Kxx,Kyy,Kxy 

transverse  permeability:  K, 


77 
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Experimental  Permeability  Prediction 


20°  (shear 


Mold  Plates 


Sheared  Fabric 


Problem:  Many  Experiments 


leaMteformation 

fafton  of  permeability  tensor 


©2003  University  of  Delaware  All  rights  reserved 


2  July  2003 


Permeability  vs.  Vf  Measurement 


♦  Detailed  measurements  of 
permeability  vs.  Vf  by 
changing  the  compaction 
load 

4  Steady  state  1-D  flow 


Experiments  were  performed  by 
Mats  Erninger 
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VARTM  Specific  Measurement : 
Distribution  Media 


Direct  (ID) 


u--..-  rz..  i — . 

5  E3  !»t  ^  j 

~  sq 

s  e-3  * 

S3  e^~  £=  a 


i 


u..^,  i.,-~  r^q 


Distribution 

Media 


Structural 

layer 


Simacek  ONR  Workshop  -  22 


©2003  University  of  Delaware  All  rights  reserved 


2  July  2003 


VARTM  Specific  Measurement : 
Transverse  Permeability 


Outer  ellipse  fit 


Inner  ellipse  fit 


Top:  Bottom  view  of  point 

injection  into  dry  preform, 
with  blackened  corn  syrup 

Left:  PERMSTAT  set-up 
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Experimental  Validation  of  TekScan 
Sensor  and  Flow 


time 


Auto  RTM 


Experiment 


Simulation  selected 
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Lesson  II. :  Describe  and  Model  the 
Process 


>  Equations  Dependent  on  Measured  (or  Predicted)  Data 
4-  Darcy’s  Equation 


© 


V  =  — 


V p 


4-  Continuity  Equation 

>  Boundary  Conditions  (More  Data) 
4-  Pressure  at  Inlet 


>  Solution 

4-  Analytic  (Closed  Form...) 

4-  Resorting  to  Computers 

>  Post-processing 

4-  Trying  to  Make  Sense  of  the  Results 
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LIMS  GUI  and  Using  LIMS 


Part  Design/Meshing 
CAD  Software 


Preparing  Data  for  Simulation 
LIMS  Ul  (User  Interface) 


Permeability 

Measurement 


Computer  Duties: 

•  Sensing 

-  Flow  Rate 

-  Pressure 

•  Control 

-  Flow  Rate 


Permeability  and  Preform 
Properties  Computation 
PERM,  DRAPE,  ... 


•Converting  Input  Data 

•Setting  Material 
Properties 

•Creating  Gates  and  Vents 

•Creating  and  Editing 
Control  LBASICFiles 


Running 
Simulation 
LIMS  or  LIMSSLV 


•Simulating  Filling 
Procedure 

•Executing  Control 
Operations 

•Writing  Desired 
Output 


Post-processing  Results 
LIMS  Ul,  TECPLOT,  ... 
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Simulation  Capabilities  I 


LIMS 


Isothermal 


Cure  ONLY 


Temperature 

& 

CURE 


Isothermal 


Cure  ONLY 
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Simulation  Capabilities  II 


Preform  Deformation 
(Draping) 


Scripting  Control 


VT 

'■1.570 

0J553 

$527 

0-502 

0.^77 

■J~DI 

0^20 

0^01 

0.375 

0.355 


HFH  Pi'OT.Pdur €r  l  g  FILL  IS  Cfl IIk.1  tp  nn  np^nad  P11h 
REh  iii  Lh  rmm  type  Bates  or  hw.^d  at  nodes 
91  as 

PLH  Ird  and  further  argurcents 
REH  KQiF] L L  h J iQ, ha^P p  1 J ,  1  £  In 


RROC  kSFUL 
DGFEhTT  MoHmks,! 

KfML  Hl»P,Mfl«0 
LET  HayD=ARuLHENTUJ 

l£1  ho-P=AfcBLn£rmzj 
let  fjaModii=.iy!5c-a 
Drn  Hd(noHode5j 

FDR  I  I  ID  NuH^h:? 

LET  Hd  ( 1  S^AfiaJWEf4T(  i  ■+£  S 
5ETtfATT£  IHtJ(  I)  .? 

HEKT  1 
DO 

SOLVE 

rofi  1=3  TO  H^ndes 

CALL  CAM BATE  Hd[1^ 
IOT  I 

LDOP  h/HM  F  rODrJUMBFRFHPT 
ERASE  HcflddiS. !  .na«P.Hfl)< 
EHPPRK 
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Simulation  Capabilities  III 


Non-lsotropic  Problems 


Graphical  User  Interface 


— 


Midplane  temperature 


Flowfronts 


=Mai 


Eik  Epp<?n 


<*  B  u|  7  *L  J  I.i  4,  lea  E  Cs  *\f)  «  I  |SSnSu  3  %  I  ■><*- 
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Simulation  Capabilities  Examples 


Saturation  of  Tows 

Saturation 
■  1.000000 

■  0.900000 

■  0.800000 
b  0.700000 

b  0.600000 

S  0.500000  II 

—  0.400000 

—  0.300000 

b  0.200000 
b  0.100000 
1  0.000000 

FASTRAC 


Sequential  VARTM  Injection 
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Distribution  Media  Modeling 
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Distribution  Media  Modeling  I 


Distribution  Media  (2D) 


Preform  (3D) 


7^ 

7^ 

P 
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Distribution  Media  Modeling  II 
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Lesson  III. :  Sensing  and  Control 


>  Is  What  Was  Predicted  Really  Happening? 

4-  Not  Sure!  There  May  Be  Various  Disturbances 

>  Let’s  Check 

4-  Common  Sense  Feedback 
4-  Sensors  and  Sensing 
4  Flowfront  Position 
4  Pressure  and  Flowrate 
4  Temperature 


>  Never  Test  For  an  Error  Condition  You  Cannot  Handle:  The 
Control 

4-  Multiple  Injection  Gates  and  Vents 
4-  Taking  Control  Action 
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Disturbances 


Where  Do  the  Disturbances  Come  From  ? 

Preform  Deformation 
Inaccurate  Preform  Cutting 
Inaccurate  Material  Properties 

How  Do  We  Detect  Disturbances  ? 

When  Process  Fails 

When  Part  does  not  Pass  Quality  Control 
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Disturbances 


How  To  Predict  Possible  Disturbances  ? 


Common  Sense 


Experience 


Rules  of  Thumb 
♦  Very  Rare 
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Sensors:  How  to  Detect  Disturbances 


>Pointwise 
>  Lineal 

♦  Electrical 

♦  Optical 
> Visual  Input 

♦  Smartweave™ 

♦  Tekscan™ 

♦  Visual 
(Camera) 
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Control 


>Possible: 

4-  Open/Close  Gates 
4-  Close  Vents 

4-  Control  Pressure  or  Flow- 
Rate  in  Existing  Gates 

4-  Tinker  with  Injected  Resin 

>Barely  Possible: 

4-  Control  Distribution  Media 

4-  Control  Compaction 

>Not  Possible 

4-  Tinker  with  Preform  and 
Distribution  Media 
Locations 
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Can  We  Model  the  Disturbances  and 
Sensors? 


>  Variations  in  Resin  Viscosity? 

♦  Yes,  Primary  Input  Value 

>  Variations  in  Preform  Permeability  and  Porosity? 
4-  As  Above 

>  Racetracking  Channels? 

4-  Multiple  Options: 

4  High  Permeability 
4 1 D  Elements 

>  Sensor  and  Control  Modeling 

4  In  LIMS  One  Can  Access  All  the  Data  During 
Simulation  and  Change  Settings  Accordingly 
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LIMS  Ul  -  Changing  Permeability 


□MS  Alone 


LIMS  Ul 


Q  File  Edit  View  lools  Options  Window  Help 

^JgJ2£j 

□  H  0  *  ^  ^  A  V  ‘a  ft  %  A  #  J0  f 

race.dmp  | 

C  Hex  ^  Dec  C  Oct  C  Bin 


I-  Inv  r  Hyp 


0.7071 06781 1 865475244008443621 0485 
Degrees  Radians  C  Grads 

r  r  Backspace  CE  C 


fnl798,  Col  97  COL  [DOS 


Hope  you  remember  element  numbers.... 
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Modeling  Racetracking  Channel  I. 


Raise  Permeability  in  Edge  Elements! 

■■TOP’’ 


94 

95 

S6 

1 

2  1 

3 

4 

/  \ 


SCALEPERM  "TOPM00  SCALEPERM  "BOTTOM"  ,  100 


\ 
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Modeling  Racetracking  Channel  II 


Add  ID  Elements! 
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1 D  Elements  for  Racetracking 
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Simple  Adaptive  Filling:  Sensor  Model 


Simacek  ONR  Workshop  -  44 


©2003  University  of  Delaware  All  rights  reserved 


2  July  2003 


Optimization 


>  So  Far  We  Just  Tried  to  Make  It 

>  The  Real  Challenge:  How  to  make  it  Better 

>  What  does  Constitute  “Better”  (Cost  Function)  ? 

4-  Faster  or  Cheaper 
4-  Stronger  or  Better-Quality 
4-  Compromise  between  those  Two 

>  What  Can  we  Modify  (Parameters)? 

4-  Part  Design  (Off-limits  for  Processing  Only) 

4-  Material 

4-  Tooling  and  Process  Parameters 
4-  Distribution  Media 

>  How  to  Relate  the  Parameters  to  the  Cost  Function 
♦  !!!  Efficiently  !!! 
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Simple  Optimization 
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Optimal  Auxiliary  Vent 


Cost  Function: 

Weighted  average  dry 
area  size  per  each 
scenario 

Nominal  case:  24. 10 


Optimization  Problem: 

Given  the  number  of  auxiliary  vents,  find  the  vent  configuration  that  will 
minimize  the  cost  function. 


Tolerances: 

A  dry  spot  which  include  3 
or  less  nodes  is  ignored. 


An  area  of  radius  3  nodes 
around  a  vent  is 
considered  safe 


V 

\ 

y 

\ 

\ 

l 

, 

\ 

i 

. 

/ 

\ 

/ 

s 

r 
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Optimal  Auxiliary  Vent  Locations 


Exhaustive  Search 

•  551  nodes  in  the  candidate  pool 

•  1 51 ,525  possibilities  for  2  auxiliary 
vents 


Exhaustive  Search  II: 

•  1 04  nodes  in  the  candidate  pool 

•  5,356  possibilities  for  2  auxiliary  vents 


Exhaustive  Search  III: 

•  7  nodes  in  the  candidate  pool 

•  21  possibilities  for  2  auxiliary  vents 
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Results,  2-Vent  Configurations 


Exhaustive  Search 


Exhaustive  Search 


5,356  vent  configurations 


21  vent  configurations 


Best  Aux.  Vent  Configuration 
Vents  @  nodes  [1  721 

Cost  function:  16.4 


Best  Aux.  Vent  Configuration 
[68  1  ]  numbered  nodes 

Cost  function:  17.4 


Accuracy  loss:  6.1  % 

Computational 
Savings:  99.6  % 
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Optimization  of  Fill  Time:  The  Quest 
for  Effectiveness 


Triagram  GAs  Exhaustive 

S. 
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Optimization:  “Distribution”  of 
Distribution  Media 
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Can  We  Expect  Any  Problems? 
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Can  We  Help  It? 
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Conclusions  I. 


>  Filling  Process  Can  Be  Accurately  Simulated 

>  Tools  Are  Readily  Available 

>  Simulation  Allows  Design  of  Process  without  Expensive 
Trial-and-Error  Stage. 

>  Process  Control  and  Optimization  Require  Flexible 
Simulation  Tools 

4-  Repetitive  Virtual  Processes  are  Required 

>  It  Is  Possible  to  Deal  with  Process  Disturbances  via 
Automated  Control 

4-  Simulation  is  Necessary  for  Design  of  Non-trivial  Process 
Control 

>  It  is  Possible  to  Optimize  Process  Design  With  Respect  to 
Various  Cost  Functions 

4-  Simulation  is  Necessary  for  Non-trivial  Cases 
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Conclusions  II. 


>  Simulation  Tools  Depend  on  Experimental 
Methods  (or  Additional  Simulation)  for  Input  Data 

4-  Permeability  Components 
4-  Also  Resin  Data,  Porosity,  etc. 

>  Simulation,  Optimization  and  Control  Require 
Occasional  Experimental  Verification 

>  User  Expertise  Is  Still  Important 
♦  Building  Part  Model(s) 

4-  Prediction  of  Possible  Disturbances 
4-  Dealing  with  Cost  Functions 
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